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MATHEMATICAL  BACKGROUND 


USE  OF  THE  DECIBEL  SCALE 

In  this  report  gains  and  losses  will  be  expressed  In  decibels.  In 
conformance  with  common  practice  In  the  field  of  communications.  The  "bel" 

Is  named  after  Alexander  Graham  Bell,  who  Invented  the  telephone.  It  Is 
defined  as  the  logarithm  to  base  10  of  an  energy  ratio.  Since  the  bel  Is 
too  large  a unit  for  greatest  convenience,  the  use  of  the  decibel  Is 
accepted  even  though  values  may  get  to  be  over  100  decibels.  The  definition, 
then , Is 


E 

decibels  (db)  = 10  log  (1.1) 

\ 

where  frequently  E Is  output  energy  and  E Is  Input  energy. 

2 1 

Considerable  liberty  has  been  taken  with  the  use  of  the  decibel  scale 
without  much  regard  for  the  strict  definition.  In  this  report  power  and 
Intensity  will  be  used  more  than  energy,  and  there  will  be  no  hesitancy  to 
express  a power  ratio  or  an  Intensity  ratio  In  decibels  by  taking  the 
logarithm  and  multiplying  by  10.  Furthermore,  an  absolute  power  or 
Intensity  level  will  be  expressed  In  decibels  with  respect  to  an  accepted 
reference  level.  The  reference  level  for  power  Is  accepted  as  on.'  watt. 

That  of  Intensity  Is  controversial.  Since  standardization  Is  pending  at 
this  time,  a reference  level  for  Intensity  acceptable  to  most  of  our  readers 
will  be  chosen,  viz,  the  power  associated  with  a pressure  level  of  one 
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(level  of  one)  mlcrobar  when  the  pressure  is  operating  into  a purely 
resistive  load. 

The  practice  of  expressing  voltage  or  current  gains  or  losses  in 
decibels  has  sometimes  been  done  incorrectly.  We  shall  consider  voltage 


gain  in  an  amplifier.  Suppose  the  input  voltage  is  e^  across  a resistance 


R.  and  the  output  voltage  is  e across  a resistance  R . Then  power  input 


is  power  output  is  e^*/R  so  that  gain  in  db  is  correctly 


e 


expressed  as  10  log  = 10  log  — — Now,  only  if  R = R can  the  R's 

^ Ci^Ro  ^ ° e^ 

be  cancelled,  and  only  for  this  case  the  expression  becomes  20  log  — . On 

®i 


the  other  hand,  in  this  example  it  is  true  that  the  general  correct 


®o  *^i 

expression  may  be  written  20  log  — + 10  log  — . While  neither  term  by 

0 , R 

i o 


itself  conforms  in  general  to  the  definition  of  decibels,  the  two 
together  do  conform,  and  no  harm  is  done  by  referring  to  each  term  as 
"decibels"  as  long  as  both  terms  are  present.  However,  we  must  be  sure 
that  both  are  present  unless  the  second  one  is  0. 

An  example  in  which  considerable  liberty  is  taken  without  in  any  way 


invalidating  the  numerical  results  is  in  some  terms  of  the  echo  ranging 

R 


equation.  We  may  write  spherical  divergence  loss  as  20  log  — in  which 

^i 


R^  and  R^  are  output  and  input  ranges  from  a point  source.  This  expression 
is  equivalent  to  10  log  ratio  of  powers  at  the  two  ranges  only  if  divergence 
loss  is  total  loss.  Therefore,  due  consideration  must  be  given  to  other 
losses  to  assure  a numerically  correct  result.  This  is  regularly  done  by 
introducing  the  remaining  loss  terms,  also  in  "decibels." 

There  is  one  common  situation  in  which  the  values  of  the  input  and 
output  resistances  of  an  amplifier  may  be  ignored  in  computing  gain  in  db. 
This  is  when  our  Interest  is  in  gain  in  signal-to-nolse  ratio.  Here  the 
ratio  of  signal-to-nolse  powers  at  the  input  is  equal  to  the  ratio  of 
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I 

! 


I 

1 


(the  ratio  of)  the  squares  of  the  signal  voltage  and  the  noise  voltage 

because  both  appear  across  the  same  Input  resistance.  Likewise,  at  the 

output  of  the  network,  resistance  is  not  a factor  in  the  ratio  of  signal  to 

noise.  If  the  ratio  of  slgnal-to-nolse  voltages  at  the  Input  were  and 

Y 

at  the  output  were  y^,  then  gain  is  correctly  expressed  as  20  log  — • Let 
It  be  consciously  recognized  that  this  gain  is  in  slgnal-to-nolse  ratio. 
OPTIMUM  VALUES 

In  underwater  sound  theory  it  Is  frequently  desirable  to  optimize  some 
function  by  the  choice  of  the  value  of  a controllable  variable.  This 
variable  is  often  frequency.  Every  term  in  the  function  may  be  expressible 
in  terms  of  frequency.  Then  again,  it  may  be  desirable  to  hold  some  of 
the  terms  constant  regardless  of  what  frequency  is  chosen.  An  example  of 
this  is  found  in  determining  an  echo  excess  where  there  is  a dependence 
on  directivity  index.  The  particular  condition  of  interest  may  involve  a 
directivity  index  that  is  a function  of  frequency  by  using  a transducer  of 
fixed  size,  or  the  condition  may  involve  a directivity  index  that  is  kept 
constant  by  doubling  the  transducer  dimensions  whenever  the  frequency  is 
halved.  In  either  case,  the  term  in  question  may  be  considered  as  a constant 
or  as  a function  of  the  variable  as  desired  in  order  to  obtain  results 
under  the  desired  restrictions. 

The  mathematical  procedure,  then,  is  to  differentiate  the  function  with 
respect  to  the  variable  in  order  to  obtain  the  slope  and  to  equate  this 
slope  to  0,  which  is  characteristic  of  an  optimum  point.  By  this  procedure, 
for  example,  we  may  solve  for  the  optimum  frequency  at  any  range.  In 
general,  there  is  a best  frequency  for  a given  acoustic  path  and  a given 


range. 
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(a  given  range.) 

STATISTICS 

General 

Whenever  the  processes  with  which  we  are  dealing  fluctuate,  thereby 
Introducing  an  element  of  uncertainty,  the  most  rigorous  mathematical 
treatment  Is  a statistical  treatment  In  which  we  deal  with  probabilities 
of  different  results. 

Functions  with  Discrete  Values 

A probability  of  the  occurrence  of  an  event  Is  the  ratio  of  number  of 
occurrences  to  number  of  Independent  trials  when  the  number  of  trials  is 
very  large.  It  is  thus  the  average  number  of  successes  per  trial  when  the 
number  of  trials  is  very  large.  Suppose  an  event  has  a probability  p = .4 
of  occurrence.  Out  of  ten  trials  we  would  expect,  on  the  average,  four 
occurrences.  By  "on  the  average"  is  meant  that  after  ten  trials  we  might 
take  another  set  of  ten,  and  another,  and  another,  and  we  might  average 
the  results.  We  might  get,  for  example,  on  successive  trials,  three, 
four,  two,  five,  four,  respectively.  On  the  first  set  of  ten  we  would 
Judge  the  probability  to  be  3/10  or  .3.  On  five  sets  we  would  judge  the 
probability  to  be  C3  + 4 + 2 + 5 + 4)/50  = .36  so  that  we  would  be  getting 
closer  to  the  true  value  of  ,4.  We  might  have  chanced  to  hit  the 
probability  .4  on  the  nose  on  the  first  trial,  but  we  wouldn't  be  sure  that 
we  had  hit  it.  Our  confidence  in  the  result  should  Increase  with  more 
trials. 

A fundamental  law  of  probability  relates  to  the  joint  probability  of  two  or 
more  independent  events.  This  joint  probability  is  the  product  of  the  individual 
probabilities.  This  may  be  expressed  as  P “ ' ”^n*  illustrate 


(Let  us  illustrate)  this  by  the  case  of  two  independent  events.  For 

simplicity,  we  shall  consider  the  probability  of  two  heads  in  two  tosses 

of  an  unweighted  coin.  The  following  sequences  are  all  of  the  possibilities, 

and  each  is  equally  likely:  H — H,  H — T,  T — H,  T — T.  Since  these  cases  are 

equally  likely,  each  must  have  a probability  of  1/4.  Thus  H — H has  a 

probability  of  1/4,  which  is  p p = 1/2  x 1/2  and  thus  Illustrates  the 

1 2 

general  rule.  We  must  keep  in  mind  that  when  we  deal  with  probabilities, 
we  are  dealing  with  average  results.  We  do  not  say  that  on  just  four  trials 
each  of  the  permutations  will  occur  just  once.  What  we  do  say  is  that  in 
many  sets  of  four  trials,  this  will  be  the  average  result. 

The  last  example  yields  more  information  than  we  sought.  It  yields 
the  following: 

P (two  heads)  = .25 
P (one  head,  one  tail)  = .50 
P (two  tails)  = .25 

We  may  ask  why  one  head  and  one  tail  is  twice  as  likely  as  two  heads.  It 
is  simply  that  there  is  only  one  way  of  tossing  two  heads,  but  there  are 
two  ways  of  tossing  one  head  and  one  tail,  either  H on  the  first  toss  and  T 
on  the  second  toss  of  T on  the  first  toss  and  H on  the  second.  A useful 
relationship  is  evidenced  here.  If  we  expand  (p  + q)^,  we  obtain  p^  + 2pq  + q^. 
Now  we  let  p = .5  and  q = .5  so  that  the  three  terms  in  the  expansion  become 
.25,  .5  and  .25.  If  p and  q are  respectively  the  probability  of  a head  and  a 
tall,  then  each  term  in  the  expansion  is  identifiable  with  the  probability  of 
one  of  the  three  combinations. 


Now  we  take  a more  complicated  application  of  the  binomial  expansion. 


(the  binomial  expansion.)  Suppose  a weather  prediction  has  the  probability 
3/4  of  being  correct.  What  is  the  probability  of  all  combinations  of 
success  and  failure  on  three  Independent  predictions?  We  expand 

(p  + ” p*  + 3p^q  + 3pq^  + 

Putting  p = 3/4  and  q = 1/4,  we  obtain  27/64,  27/64,  9/64,  1/64,  and  thus 

P (3  successes)®  27/64 
P (2  successes  only)  = 27/64 
P (1  success  only)  = 9/64 
P (no  successes)  = 1/64 

We  may  verify  this  result  by  the  case  method.  We  may  set  down 
sixty-four  cases  of  three  predictions  in  a row.  On  3/4  of  these  the  first 
prediction  will  be  successful  and  should  be  so  tabulated  in  the  first 
column,  and  the  remaining  1/4  are  unsuccessful  and  should  be  so  tabulated. 
Let  us  say  we  indicate  a success  by  one  and  failure  by  zero.  On  the  second 
trial  3/4  of  the  previous  successes  will  be  followed  by  another  success, 
and  also  3/4  of  the  previous  failures  will  be  followed  by  a success.  All 
remaining  cases  will  be  failures.  Continuing  in  this  way,  we  may  count  up 
the  various  results  and  the  probabilities  just  solved  for  above  will  be 
found  in  agreement. 
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(found  in  agreement.) 

Bernoulli  Distribution 

The  different  coefficients  in  the  binomial  expansion  represent  a 
number  of  ways  that  n things  may  be  taken  m at  a time.  From  what  we  have 
already  observed  and/or  proved,  we  may  write  a formula  for  m successes  out 
of  n trials  as  follows: 


(1  so  that) 


(m  < 2)  = .461. 


Now  the  probability  of  P (m  ^ 2)  •=  1 - .461  = .539. 

5 

The  preceding  example  was  for  a signal  present,  and  it  will  be  noted 

chat  P (m>2)  > p^.  It  is  now  assumed  that  the  signal  was  in  combination 

with  noise.  Suppose  now  that  we  have  noise  alone  with  the  probability  of 

exceeding  threshold  at  the  range  of  the  target  of  .01,  that  is,  p = .01. 

What  is  the  probability  of  two  out  of  five  marks  at  this  particular  range 

on  noise  alone?  The  same  steps  as  before  give  p (m  > 2)  = .001.  So  we 

s 

see  that  while  the  probability  of  detection  was  improved  by  requiring  two 
out  of  five  rather  than  one  out  of  one,  the  probability  of  false  alarm  was 
decreased  by  an  order  of  magnitude.  We  should  certainly  expect  to  gain 
something  by  using  five  times  as  much  time  (five  pings  instead  of  one),  and 
in  this  case  the  gain  is  apparent. 

It  is  now  asserted  that  when  n is  large,  the  binomial  distribution  is 


approximately  given  by 


P„(m)  = 

o/zv 


(m  - m) 
-~2^ 


(1.3) 


in  which  = npq,  m = np,  n is  number  of  independent  trials,  p is 
probability  of  success  in  any  single  trial,  m is  average  number  of  successes 
in  n trials,  and  m is  any  given  number  of  successes  for  which  the  probability 
is  desired,  o is  .just  a convenient  symbol  for  the  time  being. 


We  now  take  a simple  example.  We  shall  attempt  to  find  the  probability 


1 
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I 


I 


I 

f 

I 


L. 


(find  the  Jirobability)  of  m = m in  ten  trials  when  p = ,1.  For  this  case 

O = Ao  (.1)  (.9)  = .95,  m = np  = 10  (.1)  = 1,  the  exponential  term  is 

unity,  and  we  have  — ^ — = .42.  The  true  value  is  given  by  the  binomial 

oArr 

distribution  formula 


P (1)  = .39, 
10 


There  is  a deviation  of  about  8%  of  .42  from  .39.  The  central  part  of 

the  curve,  or  any  part,  does  better  when  n is  larger.  For  instance,  change  p 

_ _1 

in  the  above  example  to  .5  so  that  m = 5 and  a = 1.58.  Then  we  have  = .252. 

We  compare  this  with  p (5)  = .247.  The  deviation  in  this  case  is  about  2%. 

1 0 

The  two  cases  worked  out  above  are  plotted  for  all  integral  values  of 
m in  Fig  1.1  along  with  an  example  in  which  n is  small  and  the  approximation  is 
correspondingly  slightly  poorer.  The  approximate  form  is  plotted  as  discrete 
values  connected  by  a smooth  curve.  If  in  Eq.  (1.3)  we  replace  m - m by  r - y 
and  let  r vary  continuously  instead  of  taking  on  discrete  values,  Eq.  (1.3) 
becomes 

(r  - u)^ 

2o^ 

p (r)  = e . (1.4) 

oAtt 


This  is  the  Guassian  probability  density  function. 

Guassian  PDF  and  Distribution 

With  a random  variab’e  having  a continuous  amplitude  distribution, 
p (r)  dr  is  defined  as  the  probability  that  the  value  of  the  variable  lies 
between  r and  r + dr.  The  p (r)  is  called  a probability  density  function  (PDF). 
A Guassian  distribution  may  be  defined  by  its  PDF, 


11 


" 1/2 


(r  - 


_ 00  < r < “. 


(1.5) 


A r * XJ 

A graph  of  the  PDF  with  the  variable  expressed  as  r''  = — - — is  given  in 

Fig.  2. 


The  nth  moment  of  a variable  r with  a PDF  given  by  p (r)  is  defined  as 


\ = fL  • 


(1.6) 


The  zeroth  moment  is  unity  provided  that  the  coefficient  has  been 
correctly  specified.  The  value  of  unity  will  now  be  derived  for  the  PDF 
of  Eq.  (1.5)  by  direct  integration: 


o/2? 


(1.7) 


We  let 


dr  = adx. 


Substituting  Eq,  (1.8)  in  Eq.  (1.7),  we  obtain 


M = /”  e ^ dx. 

0 


Cl.  8) 


The  same  form  of  expression 
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(form  of  expression)  in  Eq.  (1.8)  would  hold  if  the  variable  were  y.  The 


value  M may  therefore  be  expressed  as 

0 


M 


4Zo  iZo^^  ^ 


(1.9) 


In  Eq.  (1.9)  the  integral  is  over  all  the  area  in  a plane.  We  may 
change  to  polar  coordinates  for  which  r*  = x®  + y*  and  dx  dy  = r d r d 0 and 
then  rewrite  Eq.  (1.9)  as 

^ 


- ^ 0 2tt  ' 0 


r dr 


(1.10) 


The  integral  with  respect  to  0 is  1.  The  integral  with  respect  to  r is 

- c ^ I “ l.and  therefore  M »=  1.  QED. 

0 0 

The  first  moment  is  the  integral  of  the  variable  weighted  by  the 
probability  density  function  and  is  therefore  the  average  value  of  the 
function.  For  the  Guasslan  distribution  this  is 


/TO 


,oo  r 


1/2 


(r  - p)‘ 


o/2it 


dr. 


(1.11) 


/ 

I 


^ r ~ y 

Letting  r = — ^ — , we  have  dr  = cr  dr  . Substituting  and  dividing 
Into  two  integrals,  we  obtain 


'2 


^2 


M,  = / 


+"  r a 


—00 


/2¥ 


dr  + 


^ dr. 


(1.12) 


I • 

h 


i i 

<■  ■ 


The  first  integral  is  that  of  an  odd  function  and  therefore  integrates  to 
zero.  The  second  integral  is  the  same  as  , already  evaluated  as  unity. 
Multiplying  by  its  coefficient  W,  we  obtain 


(U,  we  obtain) 


-1^- 


M = p. 
) 


The  second  moment  is  most  readily  computed  when  p = 0,  and  we  shall  treat 

this  case.  The  function  M is 

2 


We  now  have 


M = ^dr. 


(1.14) 


and  then 


dv  = — T e 


V «=  - e 


/udv  = uv  - /vdu 


|-H»  . 

uv  I = 0 

• ».0O 


0 f+“  2^ 


- / vdu  = -S-  e 

4 — OO  I / —00 


(1.15) 


(We  now  hpve) 

This  is  simply  times  M given  by  Eq.  (1.7)  when  p = 0.  Therefore, 

0 

for  y =•  0 we  obtain 


M -=  a* 
2 


(1.16) 


The  second  moment  about  the  mean  is  called  the  variance.  If  M is  computed 

2 

about  0 when  the  mean  is  p,  0*  = M - p*.  The  symbol  a = + is  called  the 

2 

Standard  deviation. 

The  distribution  function  of  a random  variable  (j)  (r)  is  defined  as  the 
probability  of  the  variable  not  exceeding  a specified  value  of  r.  This  is 


<|)  (r)  = p (r)  dr 


(1.17) 


and  for  the  Guassian  distribution 


♦ <r)  = /_ 


r 1 


- 1/2 


(1.18) 


In  tabulating  4>  it  is  convenient  to  introduce  a new  variable 


r - p 


(1.19) 


In  changing  the  scale  by  a factor  — , we  must  multiply  the  new  PDF  by  O since 
dr  = Odr",  obtaining 


p (pO  . 0 ^ - 1/2 

a/iv  /2if 


(dr  = adr"*,  obtaining) 

Eq.  (1.20)  Is  the  standard  form  of  the  PDF  and  may  be  Interpreted  as  Its  value 
with  zero  mean  and  unity  variance.  Another  good  Interpretation  Is  that  r Is 
number  of  O units  of  deviation  from  the  mean.  In  transformations  of  this 
kind  the  factor,  In  this  case  o.  Is  always  derivable  by  taking  the  derivative 
of  the  old  variable  with  respect  to  the  new,  in  this  case  from  Eq.  (1.19). 

Now  the  distribution  function  Is  given  by 

^ ir')  - — g - 1/2  r * (1.21) 

The  functions  p (r')  and  <|)  (r")  are  evaluated  in  oft  repeated  tables  (see 
Reference  1,  p.  167).  In  Fig.  1.2  p(r")  is  plotted,  and  in  Fig.1.3,1  - 4>  (r'’) 

Is  plotted.  The  form  In  Fig. 13  Is  most  useful  for  probability  of  exceeding 
a threshold.  To  get  <fi  (r)  for  any  given  r,  we  first  use  Eq.  (1.19)  to  get 
r",  and  then  we  use  Fig.L3  or  tables.  We  may  note  that 

1 - 4>  (-  rO  = $(r0  (1.22) 

For  fairly  large  r " the  tail  of  the  distribution  function  is  approximated  by 


1 - 4>  (r"*)  ''' 


1_  e ■ 

r " 


p (rO 


(1.23) 


For  example,  when  r = 4,  the  approximation  gives  .000033*  and  the  true  value 
is  .000032,  according  to  tables. 

Rayleigh  Distribution 

This  is  the  probability  associated  with  a vector  in  phase  space  having  a 
certain  magnitude  when  the  x and  y components  are  independent  and  both  have 
Gaussian  distributions.  It  is  applicable  to  an  envelope  detector.  It  also  turns 
out  to  be  the  kind  of  distribution  we  have  with  a small  Bt  product  (Bt  = 1)  after 


(product  (Bt  ■=  1)  after  ) 

a linear  detector.  In  contrast  to  Guasslan  distribution,  the  Rayleigh 
PDF  is  skewed  toward  high  values  and  always  is  positive. 


Its  derivation  is  from  two  Guassian  distributions  in  phase  space  representing 
the  X and  y components  of  a rotating  vector.  The  density  function  is  the 
product  of  their  density  functions,  giving 


- (x^  +.V  ) 


(1.2A) 


Changing  to  polar  coordinates,  we  obtain 


P (r.  6)  = ^ e 


r 

2? 


(1.25) 


The  r in  Eq.  (1.25)  is  the  Jacobian, 


3x  ^ 
3r  3r 
3x  ^ 
30  30 


= r. 


(1.26) 


Because  of  Independence  from  0,  we  obtain 


P - 2? 


(1.27) 


and 


P (r)  = -^  e 
o 


r 

W 


(.1.28) 


We  replace  o by  a since  this  quantity  is  no  longer  standard  deviation  and 


write 


r :» 


(deviation  and  write) 


P (r) 


r 


r 


(1.29) 


This  integrates  readily  to  give  probability  of  being  below  r,  that  is, 
the  "distribution  function,"  as  shown  below; 


^ (r)  = P (r)  dr  = - 


20^  ir 


= 1 - e 


r _ 
2ct^ 


(1.30) 


The  first  moment  is 


.00 

= p = ^ 


r 

2a^ 


dX' 


(1.31) 


This  Integral  is  that  of  Eq,  (I.IA)  with  a substituted  for  0 and  the  new 
expression  multiplied  by  /2Tt/2ot,  the  2 in  denominator  taking  care  of  the 
different  limits  of  integration.  Therefore,  the  result  is  that  of  Eq.  (1.16) 
with  the  substitution  of  a for  a and  the  multiplication  by  /2ir/2a  giving 


M = y = a = 1.26  a. 

1 y 2 


(1.32) 


I 

1 


We  take  the  second  moment  about  r = 0 


,oo  r' 

M «=  / ^ e 

2 ■'o 


r 

20^ 


dr 


(1.33) 


I 

1^ 

i'i 

I i 

t ■ 


We  let 


u = r 


r 

dv  = e dr 

a' 


20 


(We  let) 


du  *=  2rdr 


V = - e 


r 

2a‘‘  i<“ 


r+2  r 

< A •ft 


2^  , 

re  dr  = 


- 2a^e  r = 2a^ 


O*  = second  moment  about  the  mean  = (2  - ^)  = .43a^  (1.34) 


a = .656  a. 

Substituting  a from  (1.35)  Into  (1.32),  we  obtain 


(1.35) 


y = 1.92  a. 


(1.36) 


Having  determined  y and  a,  we  now  put  p (r)Csee  Eq.  (1.29)]  in  a form 
better  adapted  for  plotting.  Let 


r = Or 


(1.37) 


J = dr/dr' = a 


(1.38) 


- ‘ 

f , ~2oF~ 

p (r  ) = ^ r e 


(1.39) 


Using  the  relationship  between  a and  a In  Eq.  (1.35),  we  obtain 


p (r')  = .43  r'  e 


- .215 


(1.40) 


This  is  plotted  in  Fig.  1.4.  In  order  to  use  this  curve,  we  determine  v'  for 
the  given  r from  Eq.  (1.37)  and  then  read  p (r")  from  the  curve. 


-22- 


(from  the  curve.) 

The  function  <{>  (r)  is  modified  for  plotting  only  to  the  extent  of  using 
Eq.  (1.34)  so  that  Eq.  (1.30)  becomes 

_ -^3r^ 

1 - 4>  (r)  = e . (1.41) 


The  graph  in  Fig.L5  is  of  1 - <j)  (r)  , probability  of  exceeding  r,  for  greater 
convenience  than  (f>  (r). 

Chi  Square  Distribution 

If  a number  of  time  series,  each  having  a normal  distribution,  are  squared 

and  then  added  together,  there  is  obtained  a Chi  Square  PDF.  The  number  of 

normal  functions  involved  here,  n^umber  of  degrees  of  freedom,  affects  the 

particular  Chi  Square  distribution  obtained,  and  we  shall  therefore  use  a 

subscript  n to  denote  the  particular  distribution  involved.  When  n = 2,  the 

Rayleigh  PDF  of  r is  obtained  with  p (x  ^)  dy  ^ giving  p (r)  dy  when  X the 

2 2 

right  is  replaced  by  r.  This  has  already  been  covered  in  some  detail.  All 
Chi  Square  distributions  involve  a variable  squared  that  is  greater  than  zero. 
That  is,  they  range  from  zero  to  Infinity.  As  n is  increased  to  very  large 
values,  e.  g.  100,  the  Chi  Square  distribution  approaches  the  normal  distribution 
The  formula  for  the  Chi  Square  follows: 


P 


(1.42) 


The  r function  of  an  integer,  a,  is  given  by 


r (a  + 1)  = al. 


(1.43) 


When  n is  odd,  we  obtain  n/2  as  a fraction.  In  this  case,  F (n/2)  is 


obtained  using 


r (1/2)  = A 


(1.44) 


and  the  recursion  formula 


^ ■ 


• ‘'4 


i 

[ 

\ 


L. 


r (a  + 1)  = ar  (a) 


(1.45) 


Fig .16  is  a graph  of  the  PDF  as  a function  of  for  n = 2,  4,  and  10. 
The  distribution  function  <()  (x^)  is  give-',  by  the  following  formula: 


(X*)  = /o  P ^X^)  dx^ 


(1.46) 


This  function  is  plotted  for  n = 2,  4,  and  10  in  Fig.1.7. 

The  continuous  distributions  considered  so  far  are  all  characterized 
by  one  factor  that  is  an  exponential  and  another  factor  that  is  some  power 
of  the  variable.  Differences  show  up  in  this  power  to  which  the  variable  is 
raised. 

Student  ^ Distribution  with  jn  Degrees  of  Freedom 

The  Student  T probability  density  function  is  given  by  the  following 
formula: 


P (t^) 


r C(n  + l)/2] 

/iin  r 

n 

2 

1:  "J 

(n  + l)/2 


(1.47) 


This  distribution  ranges  from  minus  infinity  to  plus  infinity.  The  curves 
are  somewhat  similar  to  each  other  except  that  for  the  higher  n's;  they  become 
more  peaked.  All  have  maxima  at  t = 0.  For  larger  n the  Student  T distribution 


joint  Probabilities  of  Independent  Events 

We  can  imagine  a time  series  of  n random  samples,  each  sample  being  a 
1 or  a 0.  The  I's  represent  information  calling  for  a different  response 
or  perhaps  no  response.  We  let  all  samples  be  independent  of  each  other,  but  we 
let  the  probabilities  of  each  sample  being  a 1 or  a 0 be  p or  q respectively. 
What  is  the  probability  of  exactly  m I's  being  present? 

First  we  consider  the  probability  of  m specific  I's  and  the  remaining 
n - m samples,  all  O's.  For  this  case  of  joint  probability,  we  simply  multiply 
together  the  probabilities  of  every  sample  being  what  we  have  specified. 

There  are  then  m p's  and  (n  - m)  q's  multiplied  together  to  give 


P 

n 


(m) 


spec 


m n-m 
p q 


(1.A8) 


The  subscript  "spec"  means  specific  samples  provide  the  m I's. 

If  the  question  permits  any  m I's  instead  of  a specific  set,  the 
probability  is  that  for  a specific  set  of  m I's  multiplied  by  the  number 
of  specific  sets  that  are  possible.  As  an  example,  suppose  that  we  are 
given  a probability  of  2/10  that  any  sample  is  a 1.  Let  us  find  the 
probability  that  three  samples  will  contain  any  two  I's.  A specific  pair 

J 

of  I's  is  selected  as  I's  on  the  first  and  second  samples  with  a 0 on  the 
third  sample.  The  probability  of  this  is  .2  x .2  x .8  = .032.  But  I's  on 
first  and  third  samples  and  on  second  and  third  samples  are  also  counted  so 
that  we  have  three  combinations,  each  of  probability  of  .032  or  a total 
probability  3 x .032  = .096.  ■ 


The  formula  for  this  case  is 


In  which  C*'  is  the  number  of  combinations  of  n things  taken  m at  a time, 
m 

The  whole  expression  is  the  mth  term  in  the  binomial  expansion  of  (q  + p)". 

The  case  method  to  be  explained  next  will  make  the  picture  crystal  clear.  This 
method  may  be  profitably  resorted  to  in  many  statistical  problems.  We  shall 
introduce  some  variability  here  into  the  probabilities  on  successive  trials 
p = .2,  .3,  .4  respectively.  Of  course  q = .8,  .7,  .6  respectively.  Let 
us  consider  1000  cases  the  number  being  chosen  so  that  taking  some  number 
of  tenths  three  times  will  always  give  integers.  Then  the  following  flow 
diagram  tells  the  whole  story.  Lines  slanting  up  mean  I's,  and  lines  slanting 
down  mean  O's.  After  each  slanting  line  the  number  of  cases  falling  in 
that  category  is  given. 


FIRST 

SECOND 

THIRD 

CASES 

TRIAL 

TRIAL 

TRIAL 

NO  1 

ONE  1 TWO  I'S  THREE  I'S 

84 

144 

224 

336  

336  452 


24 

36 

56 

96 


188  24 


At  each  trial  the  assigned  probability  was  applied  to  all  the  cases.  Let  us 

follow  a specific  case  of  a 1 on  trials  1 and  3 and  a 0 on  trial  two. 

First  trial  (p  = .2) ; .2  x .1000  = 200  indicated  after  upward  line. 

Second  trial  (q  = .7):  .7  x 200  = 140  indicated  after  downward  line. 

Third  trial  (p  = .4):  .4  ^ 140  = 56  indicated  after  upward  line. 

The  totals  at  the  bottom  of  the  diagram  are  sums  of  all  combinations. 

For  two  I's  out  of  three  trials  we  have  188  out  of  the  1000  cases  for  a 


probability  of  .188. 


31 


Returning  to  Eq.  (1.49),  we  next  compute  P (m)  for  five  trials  (a)  with 

p “ .1  and  (b)  with  p = .5  in  order  to  show  the  gain  in  ability  to  distinguish 

between  these  cases  by  establishing  a criterion  involving  multiple  trials. 

m 0 1 2 3 4 5 

P (m)  (p  = .1)  .592  .327  .073  .0081  .0004  

P (m)  (p  = .5)  .031  .145  .313  .313  .156  .031 

Suppose  now  that  the  criterion  for  action  is  that  there  be  three  or  more 

I's  out  of  five.  For  p = .1,  P (n  ^ 3)  = .0085  and  for  p=  .5,  P (m^3)  = .5. 

The  ratio  of  the  two  values  of  P (m  ^ 3)  is  about  twenty-four  times  the  ratio 
of  the  p's.  This  is  a considerable  gain  in  distinguishing  between  these  cases 
by  using  a multiple  sample  criterion  rather  than  a single  sample  criterion. 

If  the  criterion  is  four  out  of  five,  the  gain  is  much  greater  still,  but  more 
cases  requiring  action  are  missed. 

j The  mathematics  involved  here  is  extensible  to  any  number  of  mutually 

I exclusive  answers,  for  example,  "yes,"  "no,"  or  "maybe."  If  the  probabilities 

i on  a single  trial  are  respectively  p,  q,  and  r,  and  if  there  are  m trials,  we 


I 

) 


I 

r * 
< 


(m  trlalB,  we)  wish  to  find  the  probability  of  i specific  "yeses"  as  oh 
trials  one,  six,  seven,  and  nine,  and  j specific  "noes,"  and  m - i - J 
"maybes."  Then  the  following  expression  for  applies; 


i jm-i-j 
(P  ) ■=  p q r •'< 

m spec  ^ 


(1.50) 


An  example  is  the  following:  given  p = 2/10,  q = 5/10,  and  r = 3/10, 

we  ask  what  is  the  probability  of  three  "yeses,"  one  "no,"  and  two 
"maybes,"  each  on  specified  trials  out  of  six  total  trials.  Substituting, 
we  obtain 


(P  ) = (2/10^^  X (5/10)  (3/10)*  = 3,6  x lO'""  • (1,51) 

6 spec 


We  consider  the  basis  for  Eq,  (1,50).  Take  n sets  of  six  trials  each, 
with  n very  large.  If  we  take  as  our  specific  case  trials  one  and  five 
"yes,"  trials  two,  three,  and  six  "no,"  and  trial  four  "maybe,"  then  out 
of  our  n sets  there  will  be  np  successes  on  the  average  on  the  n trials 
one.  This  action  reduces  the  number  of  cases  that  can  fulfill  our 
requirement  to  np,  the  remainder,  n (1  - p) , having  failed  already  to 
qualify.  Of  the  np  cases,  all  of  which  survived  the  first  trial,  npq  will 
survive  the  second.  Of  these  npqq  will  survive  the  third  and  ultimately 
npqqrpq  = np*  q^  r will  qualify,  on  the  average.  This  number  divided  by 
n gives  the  probability  of  obtaining  the  specific  result  required,  and 
this  checks  Eq . (1.50). 

Next,  remove  the  specification  of  which  particular  trials  must 
succeed  out  of  each  set  and  let  any  i "yeses,"  j "noes,"  and  m - 1 - j 


F 


("noes,"  and  m - 1 - j)  "maybes"  out  of  m trials  suffice.  Then  the 
preceding  result  must  be  multiplied  by  the  number  of  specific  ways  in 
which  this  result  may  be  obtained.  We  shall  indicate  this  number  by 


m 


. The  number  of  combinations  of  "yeses"  that  satisfy  the  requirement 

is  standing  for  i "yeses"  out  of  m trials.  For  each  of  these 

there  are  m - 1 trials  per  set  remaining,  and  the  chance  of  just  j "noes" 

m - i 

In  these  remaining  cases  is  so 


- m__m_m-l 
^Ij  ■ ‘^l  ‘^j 


(1.52) 


Let  us  note  that  we  could  start  with  the  "noes"  and  arrive  at 


m_  m m-j 
^ij  " j i 


These  last  two  formulas  may  be  checked  against  each  other.  We  let 


m = A,  i « 2,  and  j = 1.  The  first  formula  gives 

A X 3 X 2 X 1 

and  the  second  formula  gives. ■= ;; rr — 7—  = -r= 

® (3  X 2 X 1)  1 (2  X 1)  1 


4x3x2xl  2x1 

(2  X 1)(2  X 1)  1 X 1 

3x2x1 


= 12, 


The  final  formula  for  this  case  is 


Pm  (1.  j)  = Ci"  ‘ p"  q"  r‘"  - " - J.  (1.53) 


m-i  i j m-i-j 
P r 


FOURIER  SERIES 

We  may  consider  any  function  of  x having  only  a finite  number  of 
discontinuities  in  a finite  range  of  x in  which  the  function  is 
defined.  Consider  at  first  a specific  range  of  x,  namely  -tt  to  +tt. 
This  function  is  expressible  in  terms  of  an  infinite  sum  of  sines  and 
cosines  with  constant  coefficients,  as  follows: 


(coefficients,  as  follows:) 


f (x)  1^  “ j ^ r,  (^4  cos  ix  + b.  sin  ix)  . (1.54) 

* -IT  i = 0 1 1 


We  proceed  at  once  to  deriving  the  coefficients.  Multiply  both 
sides  by  cos  jx  and  integrate  from  -it  to  +it. 

The  integral  on  the  left  is  f (x)  cos  jx  dx. 

On  the  right  there  is  a sum  of  integrals  of  the  forms  a^  cos  ix  cos  jx  dx 
and  a^  sin  ix  cos  jx  dx.  All  of  these  may  be  shown  to  be  zero  except 
that  when  i = j and  both  trigonometric  functions  are  cosines.  This  result 
characterizes  a set  of  functions  called  normal  functions.  For  the  last 
case,  there  is  obtained 


cos^  jx  dx  = TTaj  j 0 


(1.55) 


a dx  = 2iTa  . 

0 ' -TT  0 


(1.56) 


Thus,  equating  right  and  left  integrals  and  dividing  by  tt  or  2ti,  we 
obtain 


^ f (x)  cos  jx  dx 
IT  •'  -IT  ■’ 


j 0 


(1.57) 


a = ~ f (x)  dx. 
0 2it  •’-IT  ' ' 


(1.58) 


In  like  manner,  it  can  be  shown  that 


b,  = — f (x)  sin  jx  dx. 
j TT  •’^TT  -* 


(1.59) 


(be  shown  that) 


A simple  example  follows;  we  let  f (x)  = 


-1  -ir  < X < 0 


1 0 < X < tr 


The  work  is  simplified  by  noting  the  f (x)  is  an  odd  function.  All  odd 
functions  are  expressible  in  terms  of  the  sines  only.  We  have 


b = ^ / - sin  jx  dx  + ^ sin  jx  dx 

jir^-TT  Ti''o 


^ f sin  jx  dx  = ^ |:uE£f_J2<|^ 
IT  •'  0 IT  ' j ' 0 


= 2 fl  - (-  1)^' 

^ [ J J 


bj=0 


j = 2,  4,  6,  8, (j  even) 


'’j  “ Ttj 


j =1,  3,  5,  7, (j  odd) 


Let  us  plot  the  sum  of  the  first  three  terms  to  see  how  closely  they  approximate 
the  function.  This  is  done  in  FigJ.lO.  The  following  observations  are 
pertinent.  In  the  flat  upper  and  lower  portions  of  the  curve,  the  maximum 
deviations  from  a fit  are  about  20%,  except  near  x = 0 and  near  x = ±180°, 
at  which  values  the  development  gives  zero.  These  exceptions  are  characteristic 
of  the  generality  that  at  a discontinuity  the  development  gives  the  average 
value.  Another  generality  is  that  the  development  always  describes  a 
periodic  function  with  one  period  traversed  in  the  range  of  definition, 
which  accounts  in  our  special  case  for  the  discontinuities  at  x = ±7r  where 


-J/- 


(at  X = in  where)  the  upper  and  lov;er  values  of  the  periodic  function  average 


X X 

The  period  2ti  from  to  +11  is  easily  changed  to  - — to  + -^  by  changing 

X 

the  arguments  of  the  cosines  and  sines  so  that  v;hen  ^ ~ 2 arguments 

are  in.  This  requires  that  the  arguments  be  Eq-  (1-59)  converts 

to  the  follovjing  equation: 


f (x) 


2 . X ( i2iTx  , , . i2Trx 

Jv 


(1.60) 


Let  us  go  through  the  evaluation  of  the  coefficients  as  before. 

We  should  become  very  familiar  with  this  step  as  we  may  forget  the  Pqs, 

(1.57),  (1.58),  (1.59),  and  the  corresponding  formulas  to  be  derived  for 

the  present  more  general  case.  Multiply  both  sides  by  cos  j — ^ and 

X X 

integrate  the  right  side  term  by  term  from  - y to  + All  terms  will 
integrate  to  zero  except 


'Ll 

n 


i ^ 0 


(1.61) 


dx  = X a. 


j = 0. 


(1.62) 


Equating  these  to  the  integral  on  the  left,  we  obtain 


^ 1 f 2 

X / f (x)  dx. 


Cl. 63) 


'X  vX 

This  is  the  average  value  of  f (x)  in  the  interval  ^ ^ to  + y.  For 


38- 


It 

i 


I 


L 


(to  + 2-  For)  j 0,  we  obtain 


2 ^ 

~ 2 f (x)  cos  jx  dx 


j j«  0.  (1.64) 


In  like  manner  of  derivation,  we  obtain 


I 


f (x)  sin  jx  dx  . 


(1.65) 


In  applying  the  general  Eq.  (1.60),  we  must  first  determine  whether 

the  function  is  even,  in  which  case  only  the  a's  need  be  determined.  It 

may  be  worth  mentioning  that  every  analytic  function  can  be  expressed  as 

an  even  function  plus  an  odd  function.  Indeed,  in  the  general  expansion 

of  Eq.  (1.52)  the  constant  term  and  the  cosine  terms  describe  the  even 

part  and  the.  sine  terms  describe  the  odd  part.  Some  functions  are  odd 

except  for  a constant  term.  In  which  case  only  a of  the  a’s  needs  to  be 

0 0 

evaluated,  e.g.  f (x)  = 0 1°^  f (x)  = 1 |®.  Since  a is  always  average 

0. 

value  in  the  whole  interval  of  definition,  it  is  obviously  0.5  for  the 
example . 

While  the  Fourier  series  properly  describes  a function  defined  in  an 

Interval  X,  it  is,  in  fact,  periodic  with  period  X so  that  it  repeats  as 

X goes  through  values  outside  of  the  interval  X. 

If  we  substitute  for  period  X,  the  reciprocal  of  frequency,  and  for  x, 

the  variable  t,  the  arguments  take  the  form  iZTif  t.  In  most  applications 

1 

t will  be  interpreted  as  time.  The  set  of  frequencies  if  (i  = 0,  1,  2***) 

'KC 

constitute  the  line  spectrum  of  the  periodic  function.  If  we  put 


a.  cos  (i2iTf  t)  + b,  sin  (i2Trf  t)  = /a,^  + b . ^ cos  (i2irf  t 
1 1 i 1 1 i I 


e) 
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(If  we  put)  with 


6 = tan 


- 

appears  as  the  amplitude  of  frequency  if^  and  tan  * — - is  the 
epoch  angle.  It  is  pointed  out  here  for  analogy  with  results  brought 
out  in  treating  the  Fourier  integral  that  even  functions  have  0*^  or  180° 
epoch  angle  and  odd  functions  have  + or  - 90°  angle, 

Fourier  Transform  of  ^ Transient  Function 

The  particular  function  used  here  is  chosen  because  it  will  later 
be  given  a physical  interpretation.  This  function  is  u (t) . For  the 
present,  however,  consider  u as  simply  a mathematical  function  of  the 
variable  t.  The  Fourier  transform  of  u will  be  denoted  by  U (f)  and  is 
given  by 


u (f)  - /!:>■  <t)  e'  ‘ " dt. 


(1.66) 


in  which  1 = /-  1.  Since 


e i 2tt  ft  _ 2iT  ft  - i sin  2v  ft, 


(1.67) 


it  will  be  appreciated  that  when  u is  an  even  function,  u cos  2tt  ft  is 

/•fOQ  #00 

= 21.  On  the  other  hand,  u sin  2it  ft  is  then  odd  and 

-00  J Q 

/+“ 

= 0.  Thus,  for  u even  we  write 


U (f)  = 2 / u cos  2tt  ft  dt, 
0 


(1.68) 


and  for  u odd  we  write 


U(f)  = -2i  Jo  u sin  2tt  ft  dt. 


(1.69) 


(even  we- write)  Very  frequently  u will  be  either  even  or  odd,  and  the 
transformation  using  either  Eq.  (1.68)  or  Eq.  (1.69)  becomes  only  half  as 
difficult  as  using  Eq.  (1.67).  However,  Eq.  (1.66)  may  in  some  cases  be 
simplest  without  conversion  by  Eq.  (1.67).  Note  that  the  Fourier 
Integral  Eq.  (i.66)  introduces  no  periodicity  as  did  the  Fourier  series. 

An  Inverse  transform  exists  for  transforming  U (f ) to  u (t) . This  is 


u(t)-r:u(f)  f'  d£. 


(1.70) 


Similar  comments  on  the  function  U (f)  being  even  or  odd  to  those  made  in 
the  preceding  paragraph  regarding  u (t)  apply  here  also.  We  shall 
pursue  this  Inverse  transform  no ^further  except  to  say  that  it  is  a bit 
tricky  to  reconstruct  a transient  from  its  transform  by  means  of  the 
Inverse  transform. 

Before  directing  the  reader  to  a reference  for  very  complete  treatment, 

we  may  remark  that  frequently  f is  replaced  by  U)  and  df  by  which 

renders  the  forward  and  inverse  transformations  less  symmetrical.  Discussion 

in  terms  of  rather  than  f frequently  meets  with  favor,  but  in  this  text 

we  will  favor  talking  about  f,  which  is  frequency  in  most  of  the  problems 

with  which  we  shall  later  deal.  Now,  for  the  very  complete  coverage  of 

Fourier  transforms,  consult  Reference  2. 

As  an  example  of  the  transforma;j^lon  of  a transient  function  of  t, 

"2 

let  us  take  u (t)  = sin  2i-  f t . Since  this  is  an  odd  function, 

° _ I 

Eq.  (1.69)  is  applicable.  We  write  2 
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(applicable.  We  write) 


U (f)  *=  - 2i  / sin  2n  f t sin  2tt  ft  dt 
0 c 


cos  2tt  (f  + f ) t - cos  2it  (f  - f ) 1 1 dt 

« 0 I 


fsin  Tf  (f  + f 
1 2tt  (f  + f ) 

t A 


) T sin  tr  (f  - f ) T 


2n  (f  - f 


) T1 

H 


^sinc(f  + f^)  T - sinc(f  - f^)  t]  . 


The  value  U (f)  for  this  case  is  plotted  in  Fig.Lll  for  T = -j-  for. 

0 

which  u (t)  is  a single  cycle.  The  significance  of  the  i becomes  apparent 
when  we  associate  t with  time  and  f with  frequency.  The  function  U (f)  contains 
a component  in  any  Interval  f,  which  is  an  odd  harmonic  function  of  time 
(not  expressed  in  the  frequency  expansion) . If  we  had  transformed 

cos  2ir  f t,  there  would  have  been  no  i Involved,  and  this  would  have 

0 

indicated  every  component  an  even  function  of  time.  In  other  words,  the 

i designates  phase  shift  from  the  phase  of  the  cosine  term.  This  phase 

indicator  could  be  complex  for  some  cases,  for  example,  for  the  transform 

of  sin  (2it  f t - 0)  , which  is  neither  wholly  even  nor  wholly  odd. 

0 

We  speak  of  U (f)  as  a spectral  density,  but  in  reality,  the  absolute 
value  of  its  square  U (f)U  *(f)  is  an  energy  spectral  density,  e.g.,  the 
energy  per  Hertz.  We  shall  now  solve  for  this  as  follows: 


U (f)  U * (f)  = ^ [sinc(f  + f ) T - sinc(f  - f ) tI*.  (1.71) 
4 1 0 0 J 


As  a check  on  our  statements,  we  are  able  to  integrate  over  the  frequency 


which  is  half  the  amplitude  squared  times  the  duration,  or  average  power  times 
time,  or  total  energy,  thereby  confirming  the  statement  that  U (f)  U * (f) 
is  energy  spectrum. 

CORRECTION  FUNCTIONS 

Before  addressing  ourselves  to  the  correlation  functions,  we  shall 
outline  our  objectives  in  this  treatment.  First,  the  cross  correlation 
and  autocorrelation  functions  will  be  defined.  Then  it  will  be  shown  that 
the  average  power  may  be  obtained  from  the  autocorrelation  function.  The 
•Fourier  transform  of  the  autocorrelation  function  will  then  be  shown  by  an 
example  to  be  the  power  spectrum  that  is  converted  to  energy  spectrum 
previously  derived  by  multiplying  by  duration  T. 

For  any  continuous  function  u (t),  the  autocorrelation  function  is 
defined  by 
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(is  defined  by)  For  two  continuous  functions  u (t  ) and  u (t  ) , the 

1 2 1 

cross  correlation  is  given  by 


, , V lim  1 

11/  (T)  = ^ 

* 2 <p->oo  * 


r ^ 

o 


(t)  u (t  - T)  dt 
] 2 


(1.75) 


For  transient  functions  these  equations  are  most  meaningful  if  the 

limit  is  not  taken.  The  averaging  time  should  equal  maximum  overlap  in 

time,  which  is  the  duration  of  the  shorter  of  the  two  functions.  Integration  need 

be  carried  out  for  any  T only  to  Include  all  t for  which  the  integrand  0. 

This  is  sometimes  less  than  the  averaging  time. 

Putting  T = 0 in  Eq.  (1.74)  obviously  gives  average  power  in  a one 

ohm  resistor.  We  say  (<(/)=  average  power. 

11  0 

We  next  consider  the  autocorrelation  function  of  sin  2iT  f t between 

0 

T T 

- -r  and  + -r  i 0 elsewhere.  Then  the  product  sin  2tt  f t * sin  2n  f (t  - l) 

II  0 0 

T T 

Is  0 except  when  t lies  between  - -f  T and  y for  positive  T and  between 
T T 

- -2  and  ^ negative  T.  The  autocorrelation  function  then  takes  the  form 


T 

‘ -f 

11  ^ / T 

-7-  + T 


sin  2iT  f t • sin  2ir  f (t  - t)  dt 
0 0 


for  T plus  and 


(1.76) 


T . 

'ir 


sin  2it  f t sin  2ir  f (t  - t)  dt 
0 0 


for  T negative. 

The  same  integration  is  involved  in  both  expressions,  namely 

/ sin  2iT  f t sin  2ir  f (t  - x)  dt.  The  quantity  sin  • sin  equals 
0 0 


(sin  • sin  equals)  cos  difference  A minus  cos  sum  Z divided  by  2.  The 


1 


I 

J 


A 

I ■ 

ft 


i > 
( ■ 


arguments  Z and  A are  2ir  f (2t  - t)  and  2n  f T respectively , so 

0 0 


T cos  2t:  f T - ->  sin  2it  f (2t  - t) 

2T  0 HIT  I 0 


. (1.77) 


The  following  equations  are  obtained  by  using  appropriate  limits: 

\ 


cos  2ir  f T 
0 


sin  2it  f (T  - T) 

]«-•'> J'-'  > »' 


and 


. 


1 1 


2T 


(r 

, sin  2tt  f (T  + X)  ) 

/ cos  2tt  f X 

1 + T)  - J.T  < 0 J 

(1 


We  note  that  <1/  here  is  an  even  function  since  a negative  t in  the  second 

expression  (that  for  negative  t)  gives  the  same  values  as  the  first 

expression  with  T positive.  Moreover,  always  an  even  function. 

Another  characteristic  of  the  autocorrelation  function  is  that  it 

goes  through  a maximum  for  T = 0.  This  may  not  be  apparent  in  Eq.  (1.78), 

but  when  the  slope  is  determined  by  differentiating  with  respect  to  T,  it 

Is  found  that  this  slope  is  positive  for  slightly  negative  T,  zero  for 

zero  T,  and  negative  for  slightly  positive  x,  with  two  terms  in  the  derivative 

nicely  cancelling  each  other  at  x = 0.  The  maximum  at  x = 0 is  never  exceeded 

by  any  other  maximum.  The  function  \p  foi  u = sin  2tt  f t is  plotted  versus 

1 1 0 

X In  Fig. 1,12  for  T = a single  cycle. 

0 


.78) 


i 
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(a  single  cycle.) 

The  only  comments  on  cross  correlation  that  will  be  made  here  are 
that  usually  the  two  functions  being  correlated  are  signal  only  and  signal 
plus  noise  and  that  this  is  an  effective  processing  method  when  the 
signal  is  spread  over  a medium  to  wide  bandwidth.  For  this  case  a filter 
would  have  to  accommodate  the  whole  bandwidth  and  would  admit  the  noise 
In  this  relatively  wide  band.  However,  cross  correlation  is  equivalent  to 
predetector  filtering  in  a bandwidth  equal  to  the  reciprocal  of  the  pulse 

length  regardless  of  signal  bandwidth. 

The  Fourier  transform  of  the  autocorrelation  function  will  now  be 

computed  for  the  case  in  which  the  time  series  is  u = sin  2tt  f t,  and  its 

0 

correlation  function  is  given  by  Eq.  (1.78),  with  i as  the  time  variable. 
Since  the  whole  correlation  function  is  even,  we  may  replace  e ^ ^ ^ by 

Cjs  2it  f t;  multiply  all  terms  of  the  autocorrelation  function  by  this 
quantity,  and  take  twice  the  Integral  with  respect  to  t from  zero  to  T since 
the  overlap  of  the  two  functions  disappears  at  just  T = T. 

We  let  the  Fourier  transform  be  denoted  by  g.  This  consists  of  two 
parts: 


a 


^ (cos  2tt  f t) 
T •'  0 0 


cos  2ir  f T d T 


(1.79) 


and 


1 *r 

— sr  / sin  2tt  f (T  - x)  cos  27T  f x d x (1.80) 
/TT  I 1 ' 0 0 ' 

0 

Starting  with  a,  we  express  the  product  of  the  cosines  as  the  cosine 
of  the  difference  plus  the  cosine  of  the  sum  over  2,  giving 
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I 


I 


A 


t • 

i 

5 "'4 

I i 


\ 


(over  2,  giving) 


We  integrate  by  parts 


1 T 


2 2T 


dv  *=  ^os  2ti  (f  " £^)  T + cos  2n  (f  + f^)  xj  d T 


du 


-dx 

2T 


T 

uv  =0 
0 


- / vdu 


2T  ^ 1^- 


sin  2it  (f  - f ) X sin  2tt  (f  + f ) x"! 


2tt  (f  - f ) 
0 


2tt  (f  + f ) 
0 


sin  2ti  (f  - f ) X sin  2it  (f  + f ) 


2ir  (f  - f ) 
0 


2it  (f  + f ) 
0 


d X 


a = - 


2T 


cos  2ir  (f  + f ) X cos  2ir  (f  - f ) x 


47T^  (f  + f ) 
0 


TJ + 


4tt^  (f  - f 

0 


= + 


2T 


1 - cos  2rt  (£  + f ) T 1 - cos  2”n  (f.  - f ) 1 

4n‘=  (f  - f )' 


0 __  ___  _0 

4tt^  (f  + f 


I 

4 


sin*  7T  (f  + f ) T sin*  ir  (f  - f ) T 

7T*  (f  + f )*%*  7T*  (f  - f )*'t* 

0 0 


l^sii 


7-  sin  c^f  + f ) T + sin  c^f  - f ) T 
‘*0  0 


(1.82) 


Continuing  with  the  term  b,  Eq.  (1.80),  we  express  the  product  of 

* 

the  sine  and  the  cosine  as  sine  of  the  sum  plus  sine  of  the  difference 
over  2,  giving 


4Tt 


jsin  2tt  [f  T + (f  - f ) X 

1 + sin  2it  fi 

F T - (f  + f ) xl! 

1 L 0 0 

J L 

0 ® J) 
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(over  2,  giving) 


7 — f /sin  2tt  f T + (f  - f ) T + sin  2tt  f’  T— (f  + f ) T I d T 
Airf^T-'o)  j^o  0 0 

^ I cos  2tt  [f  T + (f  - f ) t]  cos  2tt  [f  ^ T - (f  + f ) rj  |t 


4ir  f T 
0 


2tt  (f  + f ) 
0 


cos  2n  f T - cos  2n  f T 


4ir^  (f^  - f T 
0 


(1.83) 


Letting  f = X + y and  f = x - y and  substituting,  we  have 

0 

cos  2it  f T = cos  27T  T (x  + y)  = cos  2tt  x T cos  2tt  y T - sin  2tt  x T sin  2ir  y T 


- cos  27T  f T = cos  2ir  T (x  - y)  = - cos  2it  x T cos  2it  y T - sin  2ir  x T sin  2ir  y T 
0 

f + f f - f 

Adding  and  replacing  x by  y ^y  2~^*  with  the 

product  of  the  sines  only,  and  b becomes 


- sin  IT  (f  + f ) T sin  it  (f  - f ) T 


2it^  ^ - f ^)  T 
0 


^ sinc(f  + f ) T sinc(f  - f ) T. 
/ 0 0 


(1.84) 


Adding  Eq.  (1.82)  and  Eq.  (1.84),  we  obtain  the  perfect  square 


(f)  = - r rsinc(f  - f ) T - sinc(f  + f ) tI! 

4 1,  0 0 J 


(1.85) 


This  differs  from  U (f)  U * (f)  given  by  Eq.  (1.71)  in  that 


(Eq.  (1.71)  in  that) 


g (f)  = “ U (f)  U * (f)  . 


(1.86) 


I 


It  is  therefore  the  power  spectrum  of  u (t). 

By  far  the  easiest  way  to  obtain  g (f)  is  by  way  of  u (t)  and  Eq.  (1.76). 
However,  the  relationship  between  autocorrelation  function  and  power  spectrum 
is  quite  useful  for  general  manipulation. 

In  the  course  of  this  development  we  have  obtained  average  power 
by  squaring  the  time  series.  Integrating  over  its  duration,  and  dividing  by 
its  duration.  We  have  also  seen  that  the  power  spectrum  can  be  integrated 
from  to  +°“  to  give  average  power.  A statement  of  the  equality  of  these 
two  expressions  for  average  power  is  Parseval's  theorem. 

It  has  been  brought  out  that  no  other  maximum  in  an  autocorrelation 
function  is  greater  than  the  one  at  T = 0.  This  maximum  may  be  any  finite 
value.  For  many  purposes,  as  for  example  when  comparing  two  autocorrelation 
functions,  it  is  desirable  to  normalize  by  dividing  by  a quantity  that  will 
make  this  central  maximum  unity.  This  quantity  is,  of  course,  the  function 
itself  evaluated  at  T = 0. 

Cross  correlation  functions  generally  need  to  be  normalized  to  bring 

out  a part  of  their  significance.  This  is  done  by  dividing  if)  (t)  by 

12 

4 (oTl^  W with  ^ and  i|i  not  normalized.  When  this  is  done,  the 
11  22  11  22 

normalized  value  will  never  exceed  unity.  Its  maximum  value  is  a measure 
of  the  coherence  of  the  two  time  functions,  unity  meaning  perfect  coherence 
or,  apart  from  a constant  factor,  identity. 


r-' 

[• 


't 

I 

f 


I 

I 


f 


GAMMA  FUNCTIONS  . 

The  r function  is  defined  for  positive  real  V by: 

Co 

r(v)  = / e"*^  t'^^dt.  (1.87) 

0 

When  V is  an  integer  n,  successive  integrations  by  parts  reduce  Eq.  (1.87)  to 

r(n)  = (n  - 1)!  (1.88) 

and  as  with  the  factorial  we  have  a recursion  formula 

r(n  + 1)  = nF(n).  (1.89) 

Equation  (1.89)  holds  for  any  n,  real  or  complex,  and  may  be  used  to  give 
either  ascending  or  descending  arguments  (e.g.  solving  for  r(-l/3)  from  r(2/3)). 

The  Eq.  (1.87)  may  be  integrated  for  V = 1/2  by  the  substitution  t = u^ , 
yielding 

I 

r(l/2)  = /jT.  (1.90) 

Then  successive  applications  of  Eq.  (1.89)  yield  F(3/2),  F(5/2),  ...  , as  well 
as  F(-l/2),  F(-3/2),  .... 

. We  record  also: 

F(2/3)  = 1.35412 


F(l/3)  = 2.67893. 


(1.91) 
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BESSEL  FUNCTIONS 

Bessel  Functions  of  the  First  Kind 

In  problems  involving  cylindrical  coordinates  (e.g.  sound  propagation  in 
stratisfied  media),  we  encounter  Bessel  functions.  Most  useful  functions  are 
solutions  to  differential  equations.  The  Bessel  functions  may  be  taken  up  as 
solutions  of  Bessel's  equation. 

Bessel's  equation  is 


Z = 0 


(1.92) 


and  a solution  is  Z = J^(z) . By  writing  Z as  a power  series,  substituting 
in  Eq.  (1.92)  and  equating  coefficients  of  the  different  powers  of  z in 
the  equation  independently  to  0,  we  find  certain  relationships  among  the 
coefficients  of  the  power  series  for  Z that  are  satisfied  by  the  following 
formula  when  V equals  a positive  integer  n. 


Jn(z) 


00 


I 

m=0 


(-I)”*  (z/2)^” 
m!  (n  + m) ! 


(1.93) 


The  value  Jj^(z)  is  called  a Bessel  function  of  the  first  kind  of  order  n and 
argument  z.  The  values  of  J^(x)  and  (x)  are  plotted  in  Fig.1.13. 

For  negative  Integral  orders  it  can  be  shown  that 

• J_^(z)  = (-1)"  Jjj(z).  (1.94) 
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Equatlon  (1.93)  may  be  expressed  in  terms  of  T functions  as  follows: 


CO 


1 

m=0 


(-1)”  (z/2)^" 
r(m  + l)r(n  + m + 1)‘ 


(1.95) 


It  develops  that  for  any  positive  v we  may  substitute  this  V for  n in 
Eq.  (1.96)  obtaining 


J (z) 


? (-1)"  U/2) 

r(m  + l)P(v  + m + 1)' 


(1.96) 


Bessel  Functions  of  the  Second  Kind 

The  solution  J^(z)  evolves  from  the  straight  forward  procedure  just 

indicated.  The  student  may  substitute  the  solution  in  the  form  given  in 

Eq.  (1.93)  into  Eq.  (1.92)  thereby  verifying  that  Eq.  (1.93)  is  a solution 

■of  Eq.  (1.92).  However,  there  must  be  a second  independent  solution  of  a 

second  order  differential  equation.  This  will  be  given  without  derivation. 

Because  v enters  Eq.  (1.92)  as  v^,  substituting  -V  forV  gives  the 

identical  equation.  Therefore,  J_^(z)  is  a solution  of  Eq.  (1.92).  When  v 

is  an  integer,  J (z)  and  J (z)  are  not  independent  as  can  be  seen  from 
-n  n 

Eq.  (1.94).  However,  when  V is  not  an  integer,  it  turns  out  that  J ^(z) 
is  linearly  independent  of  J^(z),  and  a solution 


Z = AJ^(z)  + BJ_^(z) 


(1.97) 


constitutes  a "fundamental  system"  of  solutions. 
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It  should  also  be  apparent  that  the  combination  of  and  a 

function  of  Y (z)  which  is  a linear  combination  of  J (z)  and  J (z) 
V V -V 

constitutes  a fundamental  system  of  solutions  since 


} 

J 


i i 


aJ  (z)  + bY  (z)  = aJ  (z)  + bfcJ  (z)  + dJ  (z)] 

= (a  + be)  d^(z)  + bdJ  y(z) 

= AJ  (z)  + BJ  (z) 

V -U 


A precise  relationship  between  Y^(z)  and  the  J's  is  commonly  used,  namely 


Y (z)  = 
V 


j^(z)  cos  VTT  - j_^(z) 
sin  VTT 


(1.98) 


The  Eq.  (1.98)  defines  a Bessel  function  of  the  second  kind  and  is  called  a 
Weber  function.  (Ref.  3) 

When  V is  an  integer  n,  Y^(z)  is  indeterminate  from  Eq.  (1.98).  However, 
it  may  be  defined  as  the  limit  of  the  right  side  of  Eq.  (1.98)  as  V converges 
to  n.  It  has  been  evaluated  and  is  linearly  independent  of  J^(z);  therefore, 
we  have  the  equation 


Z = AJ^(z)  + BY^(z)  (1.99) 

which  holds  for  all  real  V.  Y (X)  is  plotted  in  Fig.  1.13.  Y (0)  = 

0 0 

By  using  kz  in  place  of  z in  Eq . (1.92)  and  taking  the  special  case  of 
V = 0,  we  obtain 


. 1 ^ . 0 

dz  z dz 


(1.100) 
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a form  which  may  develop  In  the  process  of  solving  the  wave  equation  In 

cylindrical  coordinates  with  k*  Introduced  In  the  separation  of  variables. 

Since  we  have  substituted  kz  for  z,  a solution  Is  J (kz) . A fundamental 

o 

system  of  solutions  Is 


Z = AJ^(kz)  + BY^(kz).  (1.101) 

It  sometimes  develops  that  only  certain  values  of  k permit  satisfying 
boundary  conditions.  These  permitted  values  of  k are  called  eigenvalues, 
and  the  Z obtained  for  any  eigenvalue  Is  an  eigenfunction  belonging  to  its 
eigenvalue,  which  in  turn  belongs  to  it.  Boundary  conditions  may  also 
require  that  either  A or  B be  zero. 

Bessel  Functions  of  the  Third  Kind 

There  are  two  functions  of  the  third  kind,  usually  called  Hankel 
functions,  and  designated  and  These  functions  are  given  by 


Hv^<z)  = (1.102) 

and 

H^^(z)  = J^(z)  - lY^(z).  (1.103) 

Either  Hankel  function  is  a solution  of  Eq.  (1.92),  and  the  two  together 
comprise  a fundamental  system  of  solutions.  In  some  solutions  to  the  wave 
equations  we  shall  use  this  fundamental  system  of  solutions.  Specifically, 
we  shall  encounter  a solution 


Z = EH  ‘(z)  + FH  *(z) 
o o 


(1.104) 
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where  E and  F are  constants.  The  function  H *(z)  will  be  Interpreted  as 

o 

an  incoming*  wave  and  H *as  an  outgoing*  wave  so  that  the  two  together  may 

o 

comprise  standing  waves  and  in  addition  traveling  waves  (when  E / F) . 

Since  J and  Y are  real,  the  Hankel  functions  are  complex  and  conjugate 

to  each  other.  The  interpretation  of  a complex  so]ution  is  an  amplitude 

|h^*|  or  and  a phase  lead  tan^'C+Y/J) , the  plus  sign  for  and  minus 

sign  for  The  amplitude  and  phase  are  tabulated  in  Ref.  9. 

We  should  note  that  this  complex  form  is  analogous  to  our  use  of  e— 

to  stand  for  harmonic  space  functions.  As  a matter  of  fact,  the  limiting 

forms  of  J , Y , and  H as  the  argument  increases  without  limit  are  cosines, 
o o o 

sines,  and  exponentials  as  follows: 


= (2/irz)^^  cos  (z  - J), 


(1.105) 


Y^(z)  = (2/Trz)«^2  (z  - J), 


(1.106) 


H '’2(z)  = (2/TTz)'/^  e±^z  - J), 

Z^^  O H 


(1.107) 


The  phase  for  each  of  these  functions  approaches  (z  - ^)  at  large  z.  The  major 
difference  from  the  harmonic  functions  is  the  factor  , which  will  be 

associated  in  our  applications  with  divergence  loss  characteristic  of  cylindrical 
spreading. 

* The  respective  roles  of  and  are  of t^n^^eversed  arbitrarily  to  match 
an  arbitrary  choice  of  the  negative  sign  solutions  of  the  wave 

equation  which  will  follow.  We  shall  use  e and  this  convention  determines 
the  specified  interpretations  of  and 


[ - ■ 
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Modified  Bessel  Functions* 

In  physical  problems  we  sometimes  encounter  negative  and  complex  wave 
numbers.  As  an  example,  if  we  have  the  complex  wave  number 

k - k + Ik  (1.108) 

o 1 


it  follows  that 


Ikx  Ikox  -kjx 
e = e ^ e * 


(1.109) 


”*lc  1 z 

The  last  factor  is  an  exponential  damping  factor.  A factor  e in  the 
Hankel  functions  of  zero  order  with  large  arguments  would  emerge  from 
Eq.  (1.107)  if  z were  replaced  by  kz,  k given  by  Eq.  (1.108). 

Imaginary  wave  numbers  will  often  occur  in  problems  of  stratified 
media.  These  may  appear  in  the  arguments  of  a Bessel  function  in  such  a way 
as  to  render  the  argument  imaginary.  Bessel  functions  of  the  three  kinds 
are  all  given  special  designations  for  imaginary  arguments  and  are  all  called 
"modified"  Bessel  functions.  We  have  the  relationships 

I^(z)  = Jy(iz)  (1.110) 

K^(z)  = Y^(iz)  (1.111) 

L^‘’*(z)  = H^*’2(iz).  (1.112) 


It  should  be  noted  that  Imaginary  arguments  lead  to  modified  functions 
which  are  real.  Furthermore,  the  variable  z here  is  not  necessarily  depth. 

* The  meaning  of  this  terminology  should  not  be  confused  with  the  meaning  of 
"Modified  Hankel  functions"  to  be  Introduced  later  on  Page  62. 
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The  modified  Bessel  functions  which  we  shall  need  have  some  function  of 
depth  as  argument. 

Up  to  this  point  we  have  shown  various  forms  of  solutions  to  Bessel's 
equation  and  have  stated  that  if  z Is  replaced  by  kz,  being  introduced 
as  a separation  constant  of  undetermined  value,  it  may  develop  that  k can 
take  on  only  certain  permitted  values  (eigenvalues)  and  still  permit 
satisfying  boundary  conditions.  Even  so,  there  is  a solution  (eigenfunction) 
for  each  eigenvalue,  and  a linear  combination  of  eigenfunctions  is  a 
solution.  Thus,  Eq.  (1.104)  as  an  example  might  more  generally  take  the 
form 

00 

Z=  I E H^‘(k  z)  + F.H  *(k.z).  (1.113) 

1^1  i o i i o i 

If  there  are  no  restrictions  on  k,  the  summation  is  replaced  by  an 
integration  over  k. 

Solutions  of  Equation  for  Depth  Function 

The  following  equation  occurs  in  the  study  of  sound  propagation  in 
stratified  media 


^ - “J'f  = 0-  (1.114) 

There  are  at  least  two  approaches  to  its  solution.  The  first  approach  which 
we  shall  describe  reduces  Eq.  (1.114)  to  a Bessel's  equation  of  order  (1/3). 
Eq.  (1.114)  is  of  the  form 


+ 


(Bz  + C)(J) 


0. 


(1.115) 


liiiiiuiuttiattiiiUK 
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We  let 


Bz  + C = (|b^)  2/3 


(1. 


and 


4.(2)  = (Bz  + C)*^2y(Q  = (|b^)3/3v(?).  (1. 


Then 


C = + C)2/2, 


(1 


and 


] 

i 


1 


I 


It 


We  then  find 


. (Bz  + C) 


(1. 


and 


d2({) 


dz 


(1 


Substituting  Eq.  (1,116),  Eq.  (1.117,  and  Eq.  (1.120)  worked  out,  into 
Eq.  (1.115),  we  obtain 


0 


(1 


116) 


117) 


118) 


119) 


.120) 


.121) 


which  is  the  desired  form. 
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Expressing  V as  Hankel  functions  we  have 


i 

1 


V = AH*(C)  + 

-f  + 


(1.122) 


i 

i 


Finally,  the  full  expression  for  <|)  in  terms  of  Hankel  functions  involving  z is 


4> 


(1.123) 


in  which 

Y =(k^(l  + 3z)  - (1.12A) 

*-o 

The  second  procedure  is  to  make  the  substitution 

K^(l  + Bz)  - 

n = (K2g)T7-J . (1.125) 

Substitution  of  Eq.  (1.125)  in  Eq.  (1.114)  yields  the  Stokes  equation 


0. 


(1.126) 


Solutions  to  this  equation  may  be  found  readily  by  expressing  <|)  in  a power 
series  of  ascending  powers,  substituting  in  Eq.  (1.126),  collecting  coefficients 
of  like  powers  and  equating  the  collected  coefficients  of  each  power  to  zero. 

The  coefficients  n “ If  2,  3,  ...,  are  all  expressible  in  terms  of  a^. 

This  gives  one  series  solution  of  the  following  form 


f(n)  = a^(i  - 

o 3. 


1*4  *7 
-^n 


9 


+ 


. . . ) . 


(1.127) 


* See  reference  6 
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The  coefficients  ^^n+l’  “ 1,2,3,..-,  are  all  expressible  in  terms 

of  a , giving  a second  solution 


/\  / 2i,.257  258io,  s 

;(n)  = a^(n  - + -7711  — + . . .) 


(1.128) 


Note  that  in  both  f(ri)  and  g(n)  the  signs  alternate  but  witli  negative 

argument  all  signs  are  alike  in  either  function.  All  other  coefficients, 

a , a are  zero. 

2 5 

The  Airy  functions,  which  have  been  tabulated,  are  related  to  f and  g 
as  follows: 


-i  -i 

Ai(-n)  = 772V  f(n)  - g(n) 


(1.129) 


Bi(-ri)  = f(n)  + g(n) 


(1.130) 


The  solution  of  Eq.  (1.126)  is  expressible  in  terms  of  the  Airy  functions 


4)  = E Al(-n)  + F Bi(-n) 


(1.131) 


Finally,  just  as  we  combined  Bessel  functions  of  the  first  and  second  kind 
to  form  Hankel  functions,  so  we  may  combine  Airy  functions  to  form  modified 
Hankel  functions 


h^  = 12'-''’  [Ai(-n)  - i Bi  (-n_)), 

h = 12’^®  [Ai(-n)  + i Bi(-n^. 


(1.132) 


Note  that  when  ti  is  positive,  the  sign  in  the  terms  for  Ai  and  Bi  ai ternate 
giving  rise  to  oscillating  functions.  When  ii  is  negative,  which  can  occur 
for  example  when  3 is  negative  and  z is  large,  Ai(-n)  is  a monat(inically 
decreasing  function  with  increasing  -n  and  Bi(-ri)  is  a monatomical ly  increasing 
function  with  increasing  -ri. 

The  modified  Hankel  functions  are  expressible  in  terms  of  z by  substituting 
Eq,  (1.125)  for  n*  The  expression  for  (}>  is  then 


(|)  = C h (n)  + D h (n)  . 
1 2 


(1.133) 


In  contrast  to  the  solution  in  regular  Hankel  functions,  Eq.  (1.123),  there  arc 


constants  as  coefficients  of  h and  h as  compared  to  y(z)  in  the  coefficients 

1 2 

of  h'  and  . This  makes  a neater  solution. 

“i"  + 

The  Airy  functions  and  their  derivatives  are  plotted  in  Fig.I.lA  and  Fig. 1, 15 
taken  from  Ref. 


LAPLACIAN 


The  wave  equation  will  be  derived  from  physical  principles  in  the 
treatment  of  sound  propagation.  In  Cartesian  coo'dinates , it  is 


yz  V = 


1 ^ 
dt 


(1.134) 


The  operator  called  the  Laplacian  needs  to  be  defined.  The  definition  will 
be  in  terms  of  Cartesian  coordinates  and  we  shall  then  show  how  to  express 
it  in  any  coordinate  system. 

The  operator  V is  defined  in  vector  analysis  as 


„ . 9 L • ^ X i.  ^ 

dx  dy  9z 


(1.135) 


In  which  1,  j,  and  k are  unit  vectors  In  the  x,  y,  and  z directions 
respectively.  The  dot  product  of  7 by  itself  is 


7^ 


(1.136) 


The  Eq.  (1.136)  defines  7^  in  Cartesian  coordinates.  When  7^  operates  on 
a scalar  quantity  4',  we  have  the  left  member  of  the  wave  equation,  Eq.  (1.13A). 
Use  of  Eq.  (1.136)  in  Eq.  (1.13A)  gives  the  wave  equation  in  Cartesian 
coordinates. 

An  expression  for  7^  in  generalized  coordinates,  q , q , and  q , is 

1 2 3 

derivable  (e.g.  see  Ref. 5)  and  takes  the  following  form: 


8 

3q 


h 

3 1 


_9. 

9q 


■)  + 


_9_ 

8q 


( 3- 

'■h  h 
1 2 


i. 

9q 


— 

-) 

»- 


(1.137) 


The  h's  are  given  by 


9q 

h,  - sr 

in  which  n is  the  normal  in  the  positive  direction  to  the  level  surface 
in  q^^,  i.e.,  the  surface  obtained  with  q^  constant. 

In  spherical  coordinates  we  have  a radial  distance  r from  the  origin, 
the  colatitude  0,  and  the  longitude  (J)  measured  from  the  positive  x axis 
toward  the  positive  y axis.  The  level  surface  in  r is  a sphere.  The 

3 IT 

normal  to  this  is  radial  in  direction  and  tt"  = h =1. 

dn  J 

The  level  surface  in  0 is  a cone,  and  the  normal  to  this  is  a ^Vector 

99  1 

normal  to  the  cone  giving  ^ = h . A level  surface  in  (J)  is  a plane  and 
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i 

1 


dn 


r sin  O' 


Using  these  h's,  Eq.  (1.137)  for  the  Laplacian 


becomes 


r^sin  0 ® 3r^  ‘‘‘  30^®^"  ® ^ 


30^ 


3(1)  sin  0 3(1) 


).  (1.139) 


For  a spherical  wave  in  which  spheres  are  equiphase  surfaces  of  constant 
amplitude,  this  reduces  to 


7^  3r  ^ 3r  ^ r 3r' 


(1.140) 


A third  system  of  coordinates  with  which  we  should  have  familiarity  is 

that  of  cylindrical  coordinates.  Here  the  three  coordinates  are  r,  the 

radial  distance  from  the  polar  axis;  z,  the  distance  from  the  origin  along 

the  polar  axis;  and  6,  the  longitude  measured  from  the  positive  x axis.  The 

3r 

level  surface  in  r is  a cylinder,  the  normal  is  radial  and  ^ = h^  = 1. 

The  level  surface  in  z is  a plane  normal  to  the  polar  axis,  the  normal  is 

3 z 

along  the  polar  (z)  axis,  and  ^ = h =1.  The  level  surface  in  0 is  a plane 

wH  2 

Including  the  polar  axis,  the  normal  is  horizontal  and  tangent  to  a cylinder 

30  1 

of  radius  r about  the  polar  axis  and  = h = — . 

3n  3 r 

Applying  the  general  Eq.  (1.138)  with  these  h values,  we  obtain 
the  Laplacian 


1 fa  9 ^ ^ 1 9^ 

r 9r  9r  9z^  r 30^ 


9_  9^1^ 

r 3r  ^ 3r  3z^  r^  30 


(1.141) 


For  cylindrical  symmetry  in  which  a sound  wave  is  traveling 
inward  or  outward  with  equiphase  surfaces  which  are  infinite  cylinders, 
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